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ABSTRACT 

The effect of gravity levels up to 400 times the local 
gravity value (a/g) was experimentally investigated in a 
copper, rotating, cylindrical boiler. Heat flux levels up 
to 80,000 BTU/hr-ft^ were reached. 

Observed improvement in heat transfer at low heat flux 
values with increasing gravity levels confirmed the results 
of previous investigations. At high heat flux values, a 
consistent trend was not observed. With increasing gravity 
levels, some data indicated a decrease in heat transfer 
while other data indicated an increase in heat transfer. 

Changes in heat transfer with variation of the liquid 
annulus thickness was noted. 
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I . 



INTRODUCTION 



The study of the effect of increased gravity levels on 
nucleate pool boiling has been underway for approximately 
eleven years. This study has been pursued because boilers 
and vapor-generation devices are being extensively revised 
because of new requirements for their use in spaceflight 
applications and in the widening use of nuclear reactors 
as energy sources. These boilers are required to operate 
at high heat flux levels and under variable gravity levels. 
The study of the effects of high gravity levels on the 
boiling phenomenon may also lead to a new understanding of 
the boiling heat transfer mechanism. 

The objective of this investigation was to obtain im- 
proved data on nucleate boiling of water at high heat flux 
levels and at gravity levels up to 400 times the local 
gravity value (a/g) . This was done in order to determine 
if there was a significant variation in heat transfer rates 
in this region with gravity level change and if there was 
any reversal in the direction of the variation. 

Costello and Tuthill [1] studied the influence of in- 
creased gravity levels of 1 to 45 a/g upon nucleate boiling 

heat transfer at heat flux levels between 100,000 and 

2 

200,000 BTU/hr ft . In this limited range of study, in- 
creased gravity levels required larger values of superheat 
for a constant heat flux. Heat transfer was decreased for 
increases in acceleration in this region. 
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Merte and Clcark [2] studied the influence of increased 

acceleration, from 1 to 21 a/g, upon nucleate boiling heat 

transfer in saturated water. Under pool boiling conditions, 

they detected an acceleration effect which was dependant 

upon the heat flux level. For values of heat flux less 

2 

than 50,000 BTU/hr ft , increased acceleration decreased 

the superheat required. At heat flux levels above 50,000 
2 

BTU/hr ft , increased acceleration increased the superheat 
required. They postulated two different modes of heat 
transfer to explain this; nucleate boiling and natural 
convection heat transfer. At low levels of heat flux, 
natural convection heat transfer governed heat removal and 
increased gravity levels enhanced heat transfer rates . At 
high heat flux levels, nucleate boiling heat transfer was 
predominant and increased gravity levels decreased heat 
transfer rates. 

In the natural convection region Ref. [3] recommends 
using formula 

Nu = C (GrPr ) m , 10 5 < (GrPr) < 2xl0 7 



for horizontal flat plates where C = .54, m = .25, Gr = 
Grashof number, Pr = Prandtl number, and Nu = Nusselt 
number. Grashof number = and is thus directly pro- 

v2 

portional to the gravity level. This shows that as gravity 
level increases, heat transfer rates in the natural convec- 
tion region also increase. 
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Adelberg [4] compared gravity forces with three other 
forces which also influence the heat transfer rate in 
nucleate boiling. These forces were bubble dynamic forces, 
surface tension, and forced convection induced drag. He 
then reviewed existing literature for experimental data 
which he considered relevant and compared his analysis 
with this data. 

Later, Adelberg [5] developed three equations for heat 
transfer in pool nucleate boiling. The range of applica- 
bility of each equation depended upon the gravity level. 
From experimental data available, he concluded that his 
equations for low and moderate a/g' levels agree favorably 
with the limited data available. No data -was available to 
verify his equation for high a/g levels. 

Graham and Hendricks [6] investigated the influence of 
acceleration upon the mechanism of heat transfer in satu- 
rated liquid boiling from 1 to 10 a/g. They demonstrated 
that increasing acceleration tended to shift the boiling 
characteristic curve in the direction of improved heat 
transfer rates in proportion to the 1/4 power of the dif- 
ference in wall temperature and fluid bulk temperature. 

Gray, Marto, and Joslyn [7] investigated heat transfer 

behavior in a rotating boiler using acceleration up to 200 

2 

a/g and heat fluxes up to 505,000 BTU/hr ft . This study 
confirmed that at low heat fluxes, increased acceleration 
increased heat transfer rates; whereas at high heat fluxes, 
the opposite was true. This was in agreement with previous 
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experimenters. The region where this crossover effect took 
place occurred at progressively higher heat flux levels for 
higher accelerations. They concluded that unresolved ef- 
fects such as surface condition and aging affected the data 
at high heat fluxes as much as did gravity variations from 
1 to 200 a/g. 

Marto and Gray [8] investigated high acceleration and 
heat flux effects on nucleate boiling. They achieved ac- 
celerations up to 400 a/g and heat fluxes to 818,000 BTU/hr 
2 

ft . Their results showed that for the inception of nucleate 
boiling from natural convection, heat flux and required wall 
superheat increased as acceleration increased. They also 
verified that at low heat fluxes, increased acceleration in- 
creased heat transfer rates. But at high heat fluxes, heat 
transfer rates were practically independent of acceleration 
variations. In this region changes in heat transfer rates 
were as small in magnitude as the usual scatter in experi- 
mental data. 

Judd and Merte [9] studied the heat flux predictions of 
six different models proposed to explain the nucleate boiling 
phenomenon. They showed that none of the models adequately 
predict nucleate boiling heat flux for other than saturated 
boiling conditions at standard gravity. All the models pre- 
dicted decreasing values of heat flux as acceleration in- 
creased. These observations suggested the omission of some 
mechanisms which are inoperative under standard conditions 
of gravity, but very sensitive to the level of subcooling 
and acceleration. 
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They presented an empirical procedure for correlating 
the results used in' their study. Their equation for heat 
flux is 



Q/A = 
+ 2.39 



0.16 Pr 1/3 (2l| 3 ) 1/3 (T wall -T B ) 4/3 (a/ g ) 1/3 



, n 6 -4.1 k(T vall T B ) 

x 10 Pr 



0 - 



9 (e t -e v ) 



(T ,,-T . ) 

wall satv 

hfg 



Metric units must be used in this equation. 

The first term is the contribution of natural convection 

and the second term the contribution of nucleate boiling. 

They showed that as a/g increases, (T^ a ^-T^) decreases when 

both (T ,,-T , ) and Q/A are constant. At low heat flux 

wall satv 7 ' 

levels, the first term is important since natural convection 

dominates. As a/g increases, the heat flux may either in- 

1/3 

crease or decrease depending upon the magnitude of (a/g) 
compared to (T , , -T ) 4/ ^ 3 . At high heat flux levels, the 
first term becomes insignificant and the second term domi- 
mates. In this region Q/A will decrease with increased a/g 
since the (T . . -T„ ) term decreases. This suggests that 
there may be a crossover region of heat flux where the ef- 
fect of increasing gravity level is reversed. 
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II. EQUIPMENT CONFIGURATION 



To conduct this high gravity nucleate boiling investi- 
gation, the boiler section of an existing rotating thermo- 
syphon apparatus was used. 

Figure 1 is a schematic diagram of the thermosyphon 
assembly and associated equipment. Figure 2 is a cross 
sectional drawing of the thermosyphon. Figure 3 is a 
photograph of the equipment setup in the laboratory. 

The equipment was the same as that used by Woodard [10] 
with the following modifications. The entire thermosyphon 
was manufactured from pure copper stock with the exception 
of the two end plates of the evaporator. These were manu- 
factured from 304 stainless steel. The grooves were designed 
to reduce heat transfer from the barrel of the evaporator 
to the end plates. The difference in the thermal conduc- 
tivity of the two different materials served the same pur- 
pose in this design. The condenser end was left open to the 
atmosphere by removing the . 25 inch threaded plug in the end 
section . 

A . EVAPORATOR MANUFACTURE 

The evaporator section was manufactured from pure copper 
stock and the end sections from 304 stainless steel. It had 
an inside diameter of 3.125 inches and an outside diameter 
of 3.635 inches. The length of the copper section was 3.125 
inches, and when assembled the overall length was 3.75 inches. 
The actual heated length was 3.5 inches. The heating coil 
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was 11-gauge Chromel-A heater wire in a 0.1875 inch inconel 
sheath helically wound around the evaporator. 

The procedure for bonding the stainless steel end plates 
to the copper section and bonding the heater wire to the 
evaporator was an exacting process. Its purpose was to en- 
sure a uniform heat flux axially along the heater surface 
and to prevent electrical shorting between the heater wire 
and other components. The procedure was as follows: 

1. A thin sheet of gold copper foil was tack welded 
to the outside evaporator surface and each end 
section. 

2. The heater wire was wound over the gold-copper 
foil with gold-copper wire wound intermittently 
between the coils. 

3. The stainless steel end plates were then pressed 
in place. 

4. The entire assembly was brazed in a hydrogen 
atmospherically-controlled oven. 

5. Sauereisen cement was air dried around the out- 
side of the heater wire. 

6. Asbestos insulation was wrapped around the heater 
wire and tied off with glass tape. 

Power was supplied to the heater coils by a 130 volt 
D.C. source capable of delivering a maximum of 100 amperes. 
This power was transmitted by a brush assembly consisting 
of eight spring-loaded graphite brushes through two bronze 
collector rings to the heater wire. The bronze collector 
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rings were supported by twelve 0.25 inch stud bolts. 

Phenolic spacers were used to align the rings and to in- 
sulate the bolts from the rings. 

The outside face of the evaporator was closed by two 
pyrex glass view ports. The inner port rested on a 
neoprene 0-ring and was separated from the outer port by 
a compressed fiber gasket. A glass tube was fused to the 
inner port and passed through a hole in the outer port. 

This tube was used for evacuating and filling the system 
when used as a thermosyphon. 

In order to measure the temperature gradient in the 
evaporator wall, four holes, 0.0625 inches in diameter and 
2.375 inches deep, were drilled radially in the evaporator 
wall. Four sheathed thermocouples were inserted into these 
holes. The center of one hole was located 0.0625 inches 
from the inner wall and another hole's center was located 
the same distance from the outer wall. These holes were 
located 180 degrees apart. The remaining two holes were 
drilled so that their centers were 0.125 inches from the 
inner wall. These two holes were also located 180 degrees 
apart. 

Due to the relative softness of copper, these holes were 
not drilled true, and it was thus necessary to measure the 
misalignment of each hole. This was done to determine the 
exact position of each thermocouple junction. The following 
procedure was used on each hole to determine its misalignment. 
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A 0.0625 inch diameter pin 2.5 inches in length was 
hardened at one end for approximately 0 . 5 inches to pre- 
clude any bending. This pin was inserted in the drilled 
hole with the hardened end outward. Protruding from the 
hole was a 0.125 inch section. A precision square 2.375 
inches high was aligned to the edge of the drilled hole 
using the protruding pin end for alignment. When pro- 
perly aligned, the precision square was clamped in place 
and the metallic alignment pin removed. Then a Tungsten 
heliarc electrode approximately six inches in length and 
0.0625 inches in diameter was inserted into the hole. This 
electrode was chosen because it is manufactured to be pre- 
cision straight. The misalignment between the precision 
square and the tungsten electrode at the top of the block 
was measured using precision feeler gages. Since the 
electrode was assumed to be perfectly straight, the mis- 
alignment of the electrode at the top of the precision 
square was the same as at the bottom of the drilled hole. 

This procedure showed that the ends of the drilled 
holes were located as follows from the centerline of the 
evaporator: 1.610, 1.672, 1.683, and 1.770 inches, respec- 
tively. The hole nearest the inner wall was only .016 
inches from the wall. A slight protrusion was visible. 

Due to the softness of the copper, the drilling may have 
caused the wall surface to bulge very slightly. 
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B. INSTRUMENTATION 

1 . Temperature Measurements 

Temperature measurements were obtained from four, 

30 gauge, copper-constantan thermocouples. These were in- 
serted in 0.0625 inch diameter holes drilled 2.375 inches 
deep in the evaporator wall. They were sheathed in 0.0625 
inches OD stainless steel tubing. 

The manufacture of these thermocouples required a 
careful procedure. First, the outer plastic insulation 
was removed from the wires by using a low flame to melt the 
insulation. The residue was then carefully scraped from 
the inner wires using a surgical blade. Extreme care was 
used so that the insulation of the inner wires was not 
ruptured. This outer insulation was removed from a one 
foot section of each thermocouple. 

Stainless steel tubing, 0.063 inches OD, was cut 
into lengths of approximately five inches. The stripped 
thermocouple wires were inserted into the tubing and the 
ends of the wires were soft soldered to make the junction. 
With this junction made, the wires were pulled back until 
only approximately 1/4 inch remained exposed. A coating 
of Conley-Weld epoxy adhesive was placed over the exposed 
wires and junction. Then the wires were again pulled back 
until only the epoxy covered junction was exposed. The 
epoxy was used to protect the junction and to hold the 
wires securely in the tube. After the epoxy had set, fine 
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sandpaper was used to make a smooth finish and reduce the 
size of the tip to that of the tubing. 

The leads from the four thermocouples in the evapor- 
ator wall were brought down the outside of the condenser 
and through the rear condenser flange to a junction board 
at the end of the shaft. From the junction board, the leads 
v/ent inside the main shaft to a slip ring assembly which 
was attached to the main shaft by a coupling. 

The output from the slip rings was fed to a Hewlett- 
Packard 2010C Data Acquisition System. This system provided 
printed or visual output of the thermocouple readings. 

2 • Rotational Speed Measurement 

The rotational speed of the evaporator was measured 
using a magnetic pickup and a Systron Donner Counter 1013 
Series. An electronic strobe light was used as a check 

3 . Condenser Cooling Water Temperature Measurement 

The inlet and outlet cooling water temperatures to 

the condenser section were measured using a Hewlett-Packard 
2801A Quartz Thermometer. The uncertainty in these readings 
was ± . 1 °F . 

4 . Electric Power Measurement 

The power into the evaporator section was measured 
using DC voltage and amperage meters. 

5 . Cooling Water Flow Rate Measurement 

The cooling water flow rate to the condenser was 
measured using a calibrated flow rotameter. The calibration 
procedure for this rotameter is described in Ref. [11]. The 
uncertainty in these readings was ±.5%. 
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6 . Liquid Annulus Thickness Measurement 



The liquid annulus thickness was measured by a paper 
rule scribed in 1/16 inch marks attached to the viewing glass 
at the evaporator end. The uncertainty in these readings 
was ±.032 inches. 

C. OPERATING PROCEDURE 

The evaporator was filled with approximately 200 ml of 
distilled water using a syringe. Filling was accomplished 
either from the condenser end or the evaporator end thru the 
glass tube. When using the latter method, the glass tube 
was sealed after filling. Cooling water flow was then es- 
tablished to a 30% reading on the rotameter. The dirive 
motor was then energized and the desired RPM established. 

When a liquid annulus had formed, DC power was applied to 
the evaporator heater. With the desired RPM and power level 
established, data was recorded when thermal equilibrium was 
established. Waiting time between data recording ranged 
from 5 to 15 minutes after a change in variables. The fol- 
lowing data were recorded: RPM, a/g, power supplied to the 
heater, four thermocouple readings, coolant mass flow rate, 
temperature difference across the condenser, atmospheric 
pressure, and liquid annulus thickness. 

Securing the machine was accomplished by first shutting 
off power to the evaporator heater and then shutting off 
the drive motor. Lastly, the cooling water was secured. 
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Ill . 



EXPERIMENTAL RESULTS 



The original method of heat flux calculation was to use 
the four thermocouples installed in the evaporator wall to 
determine the radial temperature distribution across the 
evaporator wall. With this temperature distribution deter- 
mined, the inside wall temperature, T ^ , could be computed 
along with the heat flux, Q/A, as described in Appendix B. 

Power level was limited to 6 kilowatts. Above this 
level, steaming occurred from the condenser end and con- 
siderable vapor was lost to the atmosphere. The initial 
data run revealed that with the highest power level avail- 
able, the temperature variation across the evaporator was 
too small to be accurately determined by the thermocouples. 
This was due to the high thermal conductivity of the copper 
evaporator, the thin wall, and the limit on the power level. 
The errors in the location of the thermocouple wells de- 
scribed in section II also contributed to the uncertainty 
in the thermocouple readings. Therefore this method of 
calculation was abandoned. 

The alternate method used to calculate the evaporator 
wall temperature and the heat flux was to make a heat 
balance on the condenser section. The heat supplied to 
the evaporator equals the voltage times amperage supplied 
to the heater. The heat to the condenser section equals 
the heat supplied to evaporator minus convection and radia- 
tive heat losses to the atmosphere. The heat removed from 
the condenser section equals the heat supplied from the 
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evaporator section plus the heat generated by friction in 
the seals and bearings plus the heat generated by viscous 
dissipation in the coolant. These last two sources of heat 
generation can be combined into the heat generated due to 
rotation. An overall heat balance states that the heat 
supplied to the evaporator equals the heat removed from 
the condenser minus the heat generated due to rotation. 

The measurement of the heat due to rotation was deter- 
mined at three different RPM's. Cooling water flow was 
established at a constant rate and electric power to the 
evaporator heater was secured. When steady state 
conditions were established, the temperature difference of 
the cooling water and the mass flow rate were recorded. 

The heat removed from the condenser was calculated using 
the mass flow rate and temperature difference across the 
condenser as described in Appendix B. 

Evaporator wall temperature, wall saturation temperature, 
and heat flux were then calculated by the second method de- 
scribed in Appendix B. The specific data taken were enumer- 
ated under "Operating Procedure" . All reduced data appear 
in Appendix E. 

Data runs 1 thru 5 were taken at constant RPM with vary- 
ing values of heat flux. Data taking began at 25 G's up to 
400 G's. The liquid annulus thickness decreased from .4 
inches at run 1 to .19 inches at run 5. Data runs 6 thru 11 
were also taken at constant RPM and varying heat flux. In 
this series run G started at 400 G's and decreased to 25 G's 
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during run 11. The liquid annulus thickness decreased from 
.62 inches during run 6 to .38 inches during run 11. Data 
runs 12 thru 16 were taken at constant heat flux with vary- 
ing values of RPM. Liquid annulus thickness decreased from 
.38 inches to .25 inches during this series of runs. Data 

run 17 was taken at a constant gravity level, 25 G's, and a 

2 

constant heat flux of 80,000 BTU/hr ft . This run was taken 
to determine the effect of liquid annulus thickness of the 
heat transfer coefficient. Liquid annulus thickness started 
at .5 inches and decreased to .125 inches. Data run 18 was 
taken at 400 G's with varying values of heat flux. Liquid 
annulus thickness started at .63 inches and decreased to 
.56 inches. This run was taken to determine if there was 
any effect upon the data due to increasing the waiting time 
after each variable change to at least 15 minutes. This was 
well beyond the time that steady state appeared to have been 
established from the indicator values. 

Visual observations of the state of the liquid in the 
evaporator were made during each data run. A stroboscope 
instrument was used to light the evaporator interior. At 
low gravity levels agitated nucleate boiling took place in 
the liquid annulus at all heat flux levels. The degree of 
turbulence increased as heat flux increased. At high gravity 
levels the nucleate boiling process increased very slowly as 
heat flux levels were increased. The degree of turbulence 
was much less at high gravity levels. Any effect of liquid 
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annulus thickness variation upon nucleate boiling was not 
readily apparent from visual observations. 
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IX. DISCUSSION OF RESULTS 



Figure 4 shows the results of data runs 2 thru 5. RPM 
was held constant and heat flux was varied. Gravity levels 
started at 50 G's and progressed to 400 G's. Liquid annulus 
thickness was approximately .2 inches during this series. 
This figure shows that at low heat fluxes heat transfer is 
improved as gravity levels increase. There is no crossover 
region discernible, although the data appears to merge to- 
gether at high heat fluxes. 

Figure 5 shows the results of data runs 6 thru 9 . RPM 
was again held constant and heat flux was varied. Liquid 
annulus thickness was approximately .5 inches. Data taking 
started at 400 G's and decreased to 50 G's. Again higher 
gravity levels show improved heat transfer rates at low 
values of heat flux. In this figure there is a definite 
crossover of the 100, 200, and 400 G curves. The data at 
50 G's appears to be offset at superheats which are too 
high . 

Figure 6 shows the results of -data runs 11 thru 16. 
Liquid annulus thickness was approximately .38 inches during 
these runs. Heat flux was held constant and RPM was varied. 
There is a definite improvement in heat transfer with in- 
creasing gravity level, but no crossover region is evident. 

Figure 7 shows the effect of liquid annulus change on 
the heat transfer coefficient. As the thickness of the 
annulus decreases, the heat transfer coefficient increases. 
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This agrees with the results of the investigation of Mackenzie 
[12] and the results of Ref. 13. 

Figure 8 illustrates the effect of extending the waiting 
time at each data point from 5-7 minutes up to 15 minutes. 
This figure is a comparison between data runs 6 and 18. It 
appears that the longer the waiting time the more reliable 
the data. Even though the parameter indicators show that 
steady state has been established after approximately 5 
minutes, waiting until at least 15 minutes results in a 
larger temperature difference and a higher heat flux. Thus 
it takes much longer than is readily apparent to reach 
steady state conditions. 

Figure 9 is a comparison between data collected at 25 G's 
during these tests and data presented in Ref. 2 at 20 G's 
and in Ref. 8 at 25 G's. There is good correlation between 
the results of Ref. 2, but the data of Ref. 8 has higher 
values of superheat. 

Figure 10 illustrates the data collected at 50 G's with 
that presented in Ref. 8 at the same gravity level. The 
50 G data from Ref. 8 is offset in the direction of higher 
superheat. A possible explaination for this is different 
surface conditions and ageing effects between their test 
surface and the copper surface used here. 

Figure 11 is a comparison between the 400 G data during 
this test and the 400 G data presented in Ref. 8. There 
is good correlation between both investigations. 



23 



V. CONCLUSIONS AND RECOMMENDATIONS 



A. CONCLUSIONS 

1. During nucleate pool boiling at low heat flux 
levels, there is a definite improvement in heat 
transfer with increasing gravity levels. 

2. At higher heat flux levels some data shows a 
crossover trend while other data appear to just 
merge. No definite conclusion can be drawn from 
the available information. 

3. Increasing the waiting time after any variable 
change beyond the apparent time of steady state 
enhances the reliability of the data. 

4. Variations in the liquid annulus thickness 
cause changes in the heat transfer coefficient. 

B. RECOMMENDATIONS 

1. Increase the waiting time after any variable 
change to at least 15 minutes to increase the 
reliability of the data. 

2. Operate with the machine closed to the atmo- 
sphere so that steaming is prevented and heat 
flux levels can be raised above those obtained 
in this investigation. This will necessitate 

a device for measuring the temperature or pres- 
sure in the evaporator vapor space. 
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3. Provide fins on the outside of the condenser to 
increase the heat transfer rate to the cooling 
water . 

4 . Increase the evaporator wall thickness so that 
larger temperature differences can be realized. 

5. Install a permanent photographic installation 
so that the boiling phenomenon can be visually 
recorded . 
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APPENDIX A 



CALIBRATION OF THERMOCOUPLES 

The copper-constant thermocouples were calibrated using 
the Hewlett-Packard Data Acquisition System, 2010C, for 
millivolt readout. The following fixed point calibrations 
were made : 

a. Ice point of water (32.0 °F) 

b. Boiling point of water (212.0 °F at 14.696 psia) 

c. Melting point of tin (449.4 °F) 

d. Melting point of lead (621.1 °F) 

Calibration of the ice point of water was accomplished 
by using a Thermos bottle filled with crushed ice and dis- 
tilled water. Each thermocouple was inserted in this ice 
bath and three separate sets of data were recorded. The 
average value of these three data sets was used. 

For boiling point of water, a large test tube was filled 
with distilled water, and a loose fitting cork with a hole 
in the center was inserted. Each thermocouple was inserted 
through the hole in the cork to approximately 1/2 inch 
above the water level when violent boiling had occurred. 

Four separate sets of data were taken and the average value 
used. An uncertainty for each thermocouple was also cal- 
culated using the highest and lowest reading. The largest 
value was ±.009 millivolts. The standard temperature was 
corrected for atmospheric pressure variation using the 
correction formula given in Ref. [7]. 
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Calibration at the melting point of tin was accomplished 
by heating a pure tin sample in a large test tube. A small 
diameter glass tube with one end closed was inserted into 
the sample. Then the thermocouple to be calibrated was 
placed in the glass tube. A high range glass thermometer 
was also inserted in the sample. The temperature of the 
sample was then raised above the melting point and the heat 
withdrawn. The molten sample was allowed to cool until the 
freezing temperature was reached. Then for a discernible 
amount of time the temperature stabilized. At this point 
data was taken. This procedure was repeated six times 
inorder to insure accurate results. The average of the 
six data points was used and an uncertainty value was then 
calculated for each thermocouple. The largest uncertainty 
value was ±.004 millivolts. 

The same physical setup was used to calibrate the thermo- 
couples at the melting point of lead as was used for tin. 

As the lead was cooled in the molten state to its melting 
point, it did not stablize for any appreciable time. There- 
fore a tape printout was obtained using the minimum time 
interval possible. This data was then plotted on graph 
paper and the area of minimum slope was used as the calibra- 
tion point. The average reading in this area of minimum 
slope was determined and uncertainty value again calculated. 
The largest uncertainty value was ±.l millivolts. 
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Reference 14 recommends a third order equation of the 

form 

E = at + bt 2 +ct 3 

for interpolating thermocouple readings between 100 °C and 
approximately 300 °C. Using this equation and the calibra- 
tion temperatures, the constants were determined using the 
IBM 360 computer for each thermocouple. Then using this 
result, a table of thermocouple readings versus degrees 
Centigrade and Fahrenheit was calculated for each thermo- 
couple using the computer. 

The results of the interpolating equation were plotted 
every 50 °C inorder to compare these values with the ex- 
perimental data. At the three experimental points the un- 
certainty values were also plotted. The calculated results 
compared well with the experimental data. 
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APPENDIX B 



DATA REDUCTION 



A. CALCULATION OF EVAPORATOR WALL SATURATION PRESSURE 



liquid annulus 




Evaporator 

Centerline 



r . 
1 



////////////////////////// Wall 



r_^ = evaporator inner radius 

r Q = distance from centerline to liquid annulus 
h = liquid annulus thickness 
Therefore 




Squaring both sides of the equation 

2 2 . , , . 2 
r - r . - 2r . h + h 

a i i 

The assumptions in this calculation were that solid 
body rotation takes place and that the density of the two 
phase mixture could be approximated by the density of water 
at 212 °F. This type of rotation of a fluid moving as a 
solid about an axis is called forced-vortex motion [15]. 

The pressure field varies with radial and axial position. 
Since the machine was maintained in the horizontal position 
there was no axial variation in pressure, and since rota- 
tion was at high RPM, earth's gravity was neglected. The 
radial component of the pressure becomes: 
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dp 

Dr 



( 1 ) 



pi, t,2r 

gc 



2 

where p = density of . liquid at 212 °F, Lbm/ft 

JLi 

w = angular velocity of evaporator, rad/sec 
r = radial position, ft 

2 

gc = conversion factor, 32.17 f t-Lbm/Lbf-sec 

2 

P = pressure, Lbf/ft' 

Integrating equation (1) from r to r . yields: 

ci 1 



AP = 

Substituting 



2 

P L W 

2gc 
for r 





yields : 



) 



2 

p w 

AP = — — h (2r . - h) 
2gc 



or 



p w 

P . = P , + — — h (2r . 

satw atm 1 



h) 



where 



satw 

p 

atm 



2 

pressure at wall, Lbf/ft 

2 

pressure in vapor space, Lbf/ft 



B. CALCULATION OF EVAPORATOR WALL SATURATION TEMPERATURE 
The subroutine for the method of least squares on the 
IBM 360 computer was used to find an appropriate polynomial 
to simulate the water saturation pressure-temperature curve. 
The limits of 14.0 to 50.0 psia were used and a fourth 
order polynomial gave excellent correlation between actual 
data and a number of selected test points . 
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The polynomial used was as follows: 



T 



139.435 + 6.92642P 



168663P 



2 



satw 



satw 



satw 



+ . 0024795P 



3 



0000148866P 



4 



satw 



satw 



where 

T s a tw = saturation temperature at evaporator wall, °F. 

C. CALCULATION OF TEMPERATURE PROFILE 

The heat flow through the evaporator wall during steady 
state conditions can be reduced to a one dimensional problem 
if there is no axial or circumferential heat conduction. 
Under these assumptions the steady state conduction equation 
in cylindrical coordinates is 



Integrating equation (1) once gives 
dT 

k(T)r- 5 — = constant = B.. 
dr 1 

The thermal conductivity is assumed to vary linearly 
with temperature. Therefore, 




0 



(1) 



where r = radial coordinate 



T = temperature 



k (T) = thermal conductivity of copper. 



k(T) = ko + m(T-To) 



( 2 ) 
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Data on the conductivity of copper was obtained from 
Ref. 16 . In the temperature range of interest, 0 °F to 
260 °F, the conductivity of copper is a linear function of 
temperature. This data was plotted and ko and in, the 
slope, were determined graphically. 

ko = 231.77 BTU/hr-f t-°F 
m = -.02222 BTU/hr-f t-°F 2 
To = 0 °F 



Substituting equation (2) into (1) and then integrating 
again gives: 



[ko + m (T/2-To) ] T = B.^ + lnr 



The subroutine for the method of least squares on the 
IBM 360 computer was used to find the constants B^ and B^ 
Using these constants and value of the inner wall radius, 
T wall was ex '* :ra P 0 l a ted . 

D. CALCULATION OF HEAT FLUX, Q/A 
From Fourier's Law: 



Q = 2irL [ k (t ) r ~] 



or 



Q = 2ttL B 

The inside area of the evaporator is 
A = 2irr.L 

l 
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Then 



Q = Q = ^1 BTU 

A 2nr.L r. ' hr ft^ 

1 1 



E. CALCULATION OF HEAT FLUX, Q/A, (SECOND METHOD) 
The heat transfer rate out of the condenser is 



• BTU 

Q out = m C L (T 2 " V ( hF-> 



where m = coolant mass flow rate ( ^ bm ) 

hr 



BTU 



C T = coolant specific heat (• , „„ 

L ^ Lbm-°F 



) 



T^ = coolant inlet temperature (°F) 

T^ = coolant exit temperature (°F) 

Q computed was then corrected for the heat loss due 
to evaporator rotation. 



Then Q/A for the evaporator 



Q . corrected 
out 

Area of heated surface 



Q . corrected 
out 

2irr . L 

l 

F. CALCULATION OF EVAPORATOR WALL TEMPERATURE, T ,, 

' wall 

(SECOND METHOD) 

For steady state conditions, the heat flow equation is 



Q = 



27TkL (T - T , , ) 
c wall 



In 



[17] 
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where T = average of the two center thermocouple reading, 
o F 

r c = radial distance to thermocouples, ft. 

Then dividing by the inside area of the evaporator gives: 



Q/A = 



k (T - T ..) 
c wall 



r. In — 
1 r . 

1 



Solving for T wa 2 p yields: 



T ,, = T - 
wall c 



r. (In — ) 
1 r . 

i 



k 



(Q/A) 



The value of k was assumed constant at 226.7, BTU/hr-ft °F. 
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APPENDIX C 



UNCERTAINTY ANALYSIS 

In order to determine the reliability of the data col- 
lected and to determine the most probable error, an uncer- 
tainty analysis was performed utilizing the Kline and 
McClintock technique [19] . The analysis was performed on 
the data taken during run 1 at 1 KW. 



A. UNCERTAINTY IN THE CALCULATION OF HEAT FLUX 



Q/A = q" 



m cp AT 
A 



Wq" 



W* 2 

(T*> + 

m 



W 

(_££) 

cp 



■ W A 2 W 2 
'AT ' K h ’ 



( 1 ) 



W = uncertainty in subscripted quantity 
q" = 10,045 BTU/hr-ft 2 
cp = 1.005 BTU/Lbm °F 
A = .2386 ft 2 



m = 


30% 


AT = 


2.18 °F 


11 


± . 5% 


W AT 


i — i 
+! 
11 


W 

cp 


= ±.001 


> s 

11 


±.00001 



Substituting into equation (1) yields: 




. 0487 



4.87% 
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B. UNCERTAINTY IN TIIE CALCUALT I ON OF EVAPORATOR WALL 



TEMPERATURE 



T , = T - 

wall c 



r. In — — q" 
x r i 



W. 



Tw 



(W Tc ) 2 + [g- (In 



r c l 2 

— - — ) W . ] 

r i r i rl 



,31 _ r i 
1 

c 



+ W ) + ( 1 ln ET 9" ... , 

v k r rc' v r 

k 2 



r . In 2 

+ ^ v 



( 2 ) 



r. = .130 L ft 
x 

k = 226.7 BTU/hr-f t~°F 



w ri = ±.001 

r = .1398 ft 
c 

W Tc = ±2 

w k = 

q" = 10,045 BTU/hr-f t z 



w rc = ±.001 

W „ = ±.0487 

q 



Substituting into equation (2) yields 



W, 



Tw 



2.05 

T 225.7 

w 



= .0092 = .92% 
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C. UNCERTAINTY IN THE CALCULATION OF WALL SATURATION 
PRESSURE 



P . = P , + 5 ^—' h ( 2 r . -h ) 

satw atm 2g c 1 




Wp satw “/ (W Pa tm > 2 + ( 5iI (2r i' h)W P )2 



pwh' 



W^ 



( 2 ri-h)W w2 )‘ + (^-2 h( 2 r i -h)Kg c ) 



+ (^ W r .) 2 ( 2 r r 2 h ) Wh ) 2 



(3) 



r. = .1302 ft 
x 

W = 78.5 rad/sec 
Wg c = ±.001 

h = .0338 BTU/hr-f t 2 -°F 

W Patm = 11 -° 

w r . = ±.001 

r x 

g c = 32.17 ft/sec 2 

W = ± . 05 
P 

W, = ±.032 
h 

L =60 Lbm/ft 2 



W 2 = ±.l 



Substituting into equation (3) yields: 



Wp 



satw 3.36 



P , 2158.12 

satw 



= .0015 = .15% 
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D. UNCERTAINTY IN THE CALCULATION OF WALL SATURATION 
TEMPERATURE 



T , = 139.434 + 6.92642 P . - .168663 P ^ ' 

satw satw satw 

+ .0024795 P , 3 - .0000148866 P . 4 

satw satw 



9 T 



satw 



W. 



T satw 



9 P , Wp satw 
satw 



(4) 



Substituting into equation (4) yields 



W. 



T satw 



212 . 95 



= .015 = 1 . 5 ! 



satw 

The largest uncertainty appears in the calculation of 
the heat flux. This is due to the error in the mass flow 
rate measurement and to the error in the coolant flow tem- 
perature difference measurement. 
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APPENDIX D 
FIGURES 
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FIGURE 



thermocouple junction board 



i 
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FIGURE 2 . CROSS SECTION OF ROTATING HEAT P PE 



CAUTION 




FIGURE 3. 
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FIGURE 7. EFFECT OF ANNULUS THICKNESS ON HEAT 
TRANSFER COEFFICIENT 
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FIGURE 10. COMPARISION OF 50 G DATA AND 50 G DATA [8] 
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FIGURE 11. COMPARISON OF 400 G DATA AND 400 G DATA [8] 
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APPENDIX E 
DATA 
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Run 1, P atm 14.681, T satv 211.95 
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Run 3 ' P atm 14 ‘ 681 ' T satv 211 ' 95 
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